Stability of Periodic Soliton Equations under Short Range Perturbations 
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We consider the stability of (quasi-) periodic solutions of soliton equations under short range 
perturbations and give a complete description of the related long time asymptotics. 

So far, it is generally believed that a perturbed periodic integrable system splits asymptotically 
into a number of solitons plus a decaying radiation part, a situation similar to that observed for 
perturbations of the constant solution. We show here that this is not the case; instead the radiation 
part does not decay, but manifests itself asymptotically as a modulation of the periodic solution 
which undergoes a continuous phase transition in the isospectral class of the periodic background 
solution. 

We provide an explicit formula for this modulated solution in terms of Abelian integrals on the 
underlying hyperelliptic Riemann surface and provide numerical evidence for its validity. We use 
the Toda lattice as a model but the same methods and ideas are applicable to all soliton equations 
in one space dimension (e.g. the Korteweg-de Vries equation). 

PACS numbers: 05.45.Yv, 02.30.1k, 02.70.Hm 



I. INTRODUCTION 

One of the most important defining properties of soli- 
tons is their stability under perturbations. The classi- 
cal result going back to Zabusky and Kruskal [l3| states 
that a " short range" perturbation of the constant solu- 
tion of a soliton equation eventually splits into a number 
of stable solitons and a decaying background radiation 
component @ . So the solitons constitute the stable part 
of arbitrary short range initial conditions. It is generally 
believed, and is claimed in |9j], that this remains valid 
when the constant background solution is replaced by a 
quasi-periodic one. Our aim here is to show this is not 
the case, and that a thus far undiscovered phenomenon 
appears in the description of the long time asymptotics. 

Solitons on a (quasi-)periodic background have a long 
tradition and are used to model localized excitements on 
a phonon, lattice, or magnetic field background. Conse- 
quently, periodic solutions, as well as solitons travelling 
on a periodic background, are well understood. The first 
results were given over thirty years ago in the pioneering 
work by Kuznetsov and Mikhailov [9j, where the stabil- 
ity of solitons of the Korteweg-de Vries equation on the 
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background of the two-gap Weierstrass solution is inves- 
tigated. There the iV-soliton solution on this background 
is computed and it is shown that each soliton produces 
a phase shift. In addition, it is claimed that the asymp- 
totic state of any short range initial condition is a set of 
solitons. 

However, we will show that the asymptotic state is 
more complicated. The reason for this is related to the 
fact that the phase shifts of the solitons do not necessarily 
add up to zero. Hence there must be an additional feature 
making up for the overall phase shift. One might conjec- 
ture that solitons alone suffice to describe the asymptotic 
state at least in the case where the phase shifts add up 
to zero — which is the situation assumed in 0. How- 
ever, we will show that even if no solitons are present, 
the asymptotic state is not just the periodic background. 

Due to the lack of powerful asymptotic analysis tools 
and computer power at the time of 9] it was impossible to 
identify this contribution then. However, it seems that 
this omission was neither pointed out in the literature 
nor was a complete description of the asymptotic state 
known. 

To illustrate these facts, let us consider the doubly 
infinite Toda lattice in Flaschka's variables (see e.g. [ijj 
or El) 

b(n, t) = 2{a(n, tf - a(n - 1, tf), 
d(n, t) = a(n, t)(b(n + 1, t) — b(n, £)), 

(n, t) G Z x R, where the dot denotes differentiation 
with respect to time. We will consider a quasi-periodic 
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II. QUASIPERIODIC SOLUTIONS IN TERMS 
OF RIEMANN THETA FUNCTIONS 



We begin by recalling that quasi-periodic solutions are 
most conveniently written in terms of Riemann theta 
functions (@D- Such explicit formulae were first derived 
in the mid-70s thanks to the pioneering work of the Soviet 
school (Novikov, Its, Matveev, etc.). In the particular 
case of the Toda lattice the underlying Riemann surface 
is associated with the square root 
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FIG. 1: Numerically computed solution of the Toda lattice, 
with initial condition a period two solution perturbed at one 
point in the middle. 



where the numbers 



spectrum S = \J?_jE<2-i , E; 



are the band edges of the 
( nL „2.7, ^2.7+ 1] of the corresponding Ja- 
cobi operator (see [Tl| ) . We can picture this surface as 
the end result of cutting and pasting; we consider two 
"sheets" (copies of the complex plane), we cut them along 
the segments [E ,Ei], [E 2 , E 3 ] , . . . and paste the top side 
of the upper sheet to the bottom side of the lower sheet 
and the bottom side of the upper sheet to the top side of 
the lower sheet across every such segment. 

On this Riemann surface we have a standard basis of 
normalized holomorphic differentials 



algebro-geometric background solution (a q ,b q ) (e.g., any 
periodic solution) plus a short range (in the sense of 0) 
perturbation (a, b). The perturbed solution can be com- 
puted via the inverse scattering transform. The case 
where (a q ,b q ) is constant is classical (see again ^ or 
[T^I'l. but the more general case applicable here has only 
recently been analyzed in Q (see also Q). 

In Figures ^the numerically computed solution corre- 
sponding to the initial condition a(n, 0) = 1, b{n, 0) = 
(— l) n + 2<5o,n is shown. In each of these two pictures the 
two observed lines express the variables a(n,t) at two 
frozen times t = 90 and t — 165. In areas where the lines 
seem to be continuous this is due to the fact that we 
have plotted a huge number of particles (around 1000) 
and also due to the 2-periodicity in space. So one can 
think of the two lines as the even- and odd-numbered 
particles of the lattice. We first note the single perturba- 
tion which separates two regions of apparent periodicity 
on the left. Comparing the two pictures we note that 
this is a solitary wave travelling to the left. Our com- 
putations show that it is indeed a soliton. To the right 
of the soliton, we observe three different areas of appar- 
ent periodicity (with period two). Between these areas 
there are transitional regions which interpolate between 
the different period two regions. 

It is the purpose of this paper to give a complete ex- 
planation of this picture. We provide the details in case 
of the Toda lattice though it is clear that our arguments 
apply to other soliton equations as well. 



R 
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-dz, 1 < j < g 



(the constants Cj (k) have to be determined from the usual 
normalization with respect to a canonical homology ba- 
sis). 

Introduce the vector 

r°°+ 9 rfJ-k 

Zj{n,t) = Q 

•J Eo k—1 ^0 

-n C j +t2c j (g)-E Eo €C». 

J oo_ 

Here oo± are the points above oo on the upper/lower 
sheet and S Bo is the vector of Riemann constants. The 
numbers /ij are some arbitrary points whose images in 
the complex plane lie in the j-th interior spectral gap. 
All possible choices form the isospectral class of quasi- 
periodic Jacobi operators with the given spectral bands. 
This isospectral class is just a g dimensional torus, since 
the preimage of each gap with respect to the map that 
maps the surface to the complex plane (seen as one of 
the two sheets) consists of two parts, one on the upper 
and one the lower sheet, which form a circle. 

Then the well-known formulas for the solutions read 



a(n, t) 



2 -~a 2 - 



6(n,*) = 6 + ||ln 



{z(n + l,t))9(z{n - \,t)) 
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where a, b are the averages and 

0(z) = J2 ex P -'"i {ill- ±) + { -^ m) ) . ; s O" . 



is the Riemann theta function of our surface. The matrix 
r is the matrix of ^-periods of the normalized basis of 
holomorphic differentials Q. 



III. THE MAIN RESULT 

We will assume for simplicity that no solitons are 
present. In other words, we assume that the Jacobi op- 
erator above has no eigenvalues. The assumption is not 
crucial, since the solitons can be easily incorporated us- 
ing a Darboux transform. 

To obtain the long time asymptotics one reformulates 
the problem as a Riemann- Hilbcrt problem (RHP) on the 
underlying Riemann surface. This shows that the solu- 
tion can be read off from the two by two matrix valued 
function m that is meromorphic off the preimage of the 
spectrum E, with divisor satisfying 

(rriji) > («.*)> ( TO J2) ^ - Y^ V ^{n,t), 

k k 



j = 1,2, jump matrix given by 



J(n nf) - fl-\R(p,n,t)\ 2 -R(p,n,t)\ 

J{p,n,t)-y R(pnt) 1 y 

p G S, and normalization 

m(p, n, t) — * I, as p — > oo + . 

This is similar to the RHP applicable to the constant 
background case, with the main difference being that we 
allow poles at the points /ij(n,t), and their flip images 
t) on the other sheet. The g points fj,j(n,t) are 
uniquely defined by the Jacobi inversion problem 

E/ = E/ Cj-t2 Cj (g). 

k=l E ° k=l E ° Joo- 

This is a crucial point and related to the fact that our 
RHP is no longer formulated in the complex plane. While 
a holomorphic RHP with jump of index zero has a solu- 
tion in the complex plane, this is no longer true on a 
surface of genus g unless we admit at least g poles ( ) . 
In particular, the above RHP has no holomorphic solu- 
tions except in special cases (e.g. if there is no jump). 
In fact, a main issue in the mathematical analysis of the 
RHP (jgj) is to define an appropriate space of solutions 
which makes the problem well-posed. 

The matrix elements of the jump are of the form 

R(p, n, t) = R{p)S(p, n, t) exp(^(p)), 



where R{p) is the reflection coefficient at t = 0, Q is a 
ratio of four theta functions 



9(p,n,t) 



e(z(p,n,t))6(z(p*, 0,0)) 



6(z(p,O,O))0(z(p*,n,t)) 
(of modulus one) , and the phase 4> is given by 



<j>(p) = 2 n Q + 2- uj 



Here 



dz 



is the normalized Abelian differential of the third kind 
with poles at oo+ respectively oo_ and 
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n^oQs-Aj) 



R. 
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is the normalized Abelian differential of the second kind 
with second order poles at oo + respectively oo_. The 
constants Xj respectively Xj again have to be determined 
from the normalization. 

There are g + 1 stationary phase points Zj(n/t) (roots 
of 4>') which behave as follows: As r\ — j runs from 
—oo to +oo we start with z g (vi) moving from oo towards 
E2g+\ while the others stay in their spectral gaps until 
z g (rf) has passed the first spectral band. After this has 
happened, z g -\ can leave its gap, while z g (rj) remains 
there, traverses the next spectral band and so on. Until 
finally zo(rj) traverses the last spectral band and escapes 
to — oo. 

Factorizing the jump matrix and using the asymptotic 
analysis of oscillatory Riemann-Hilbert problems intro- 
duced in ( 0] ) , but generalized accordingly in Q , we have 
shown that for long times the perturbed Toda lattice is 
asymptotically close to the following limiting lattice de- 
fined by 

ai(j,t) 2 6{z{n,t)) 6(z(n + l,t) + 5(n,t)) 



a q (j,t) 9(z{n + l,t)) 6(z(n,t)+5(n,t)) 



X exp J [ !0g(l ~ |i?| 2 Voo_oo + I 

\2ir\Jc(n/t) j 



St(n,t) 



2ni 



log(l - \R\ 2 )Q, 

C(n/t) 



where R is the associated reflection coefficient, C(n/t) — 
S n (— oo, Zj(n/t)), and Zj{n/t) is the single stationary 
phase point lying in the spectrum, if there is such a point, 
or otherwise, one of the two stationary phase points lying 
in the same spectral gap. 

In other words, as in p| and one has to factor- 
ize the jump matrix and deform the contour according 
to the stationary phase points. The deformed RHP is 
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will be given in ||. Instead, we are here providing a 
numerical confirmation of the above result. Indeed the 
limiting lattice can be easily computed numerically. For 
the initial data defined in the introduction a and at time 
t = 90 the result is shown in Figure |2 The soliton was 
added using a Darboux transformation. 




FIG. 2: Solution of the Toda lattice (black) plus correspond- 
ing limiting lattice (gray) . No gray points visible implies that 
black points overlap gray points. 

asymptotically close to the unperturbed one by the os- 
cillatory nature of the jump. However, due to the poles 
at fJ,j(n, t) there is an additional contribution which gives 
raise to the modulated limiting lattice. We thus have here 
an extension of the so called nonlinear stationary phase 
method to problems defined on a Riemann surface. 

In summary, for any short range perturbation a(n, t) 
of a quasi-periodic solution of the Toda lattice one has 
that YljLn ail ft) ~ * ^ uniformly in n, as t — > oo. From 
this one recovers the a(n, t) and by differentiating one 
recovers the b(n,t). It thus follows that 

\a(n,t) — ai(n,t)\ + \b(n,t) — bi(n,t)\ — » 

uniformly in n, as t — » oo. 

If solitons are present we can apply appropriate Dar- 
boux transformations to add the effect of such solitons. 
What we then see asymptotically is additional travelling 
solitons on a periodic background or 0). 

We are not presenting here our detailed computations 
leading to the asymptotic identification of the perturbed 
problem described in the introduction and the modulated 
lattice described above. A complete mathematical proof 



IV. CONCLUSION 

Let g be the genus of the hyperelliptic Riemann sur- 
face associated with the unperturbed solution. We have 
shown that the n/t-plane contains g + 2 areas where the 
perturbed solution is close to a quasi-periodic solution in 
the same isospectral torus. In between there are g + 1 
regions where the perturbed solution is asymptotically 
close to a modulated lattice which undergoes a continu- 
ous phase transition (in the Jacobian variety) and which 
interpolates between these isospectral solutions. In the 
special case of the free solution (g — 0) the isospectral 
torus consists of just one point and we recover the clas- 
sical result. 

To conclude, we believe that, apart from the interest- 
ing mathematical problem of the analysis of Riemann- 
Hilbert problems on Riemann surfaces, the limiting lat- 
tice is interesting on its own. We even speculate that 
there may even be a physical or technological interest in 
such a modulated lattice and that the asymptotic prob- 
lem described here might provide a feasible way of pro- 
ducing such modulated solutions. 
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